Problèmes aux limites de Laplace en 2D-3D - Considérations générales — p 1 


Solutions analytiques d'équations aux dérivées partielles bidimensionnelles par la méthode de 
séparation des variables et le principe de superposition des solutions particulières dans la 
catégorie des équations de Laplace, Poisson et l'équation de diffusion (chaleur) 


Extension à quelques systèmes de coordonnées orthogonales 
Cartésien, Cylindrique, Sphérique, Ultrasphérique, Sphéroïdale allongé et aplati, Elliptique 
cylindrique, Parabolique cylindrique, Toroïdal, Paraboloïdal 
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Dans le vif du sujet : les équations de Laplace, Poisson et l'équation de la chaleur 
Conditions aux limites et conditions initiales 
Conditions aux limites, conditions initiales homogènes, inhomogènes 
Principe de superposition de l'équation de Laplace (et également régime stationnaire de 
l'équation de la chaleur) 
Principe de linéarité des solutions du problème aux limites stationnaires 
Séparation des variables 
Problème de Sturm-Liouville régulier (et auto-adjoint) 
Problème de Sturm-Liouville régulier, auto-adjoint avec des valeurs propres tendant vers - °° 
Problème de Sturm-Liouville singulier et auto-adjoint 
Calcul des normes des fonctions propres du problème de Sturm-Liouville 
Exemple : cas des fonctions de Legendre, des fonctions de Bessel, des fonctions sinusoïdales 
Finitude des solutions 
Singularité du système de coordonnées et régularité du gradient 


Solution de l'équation de Laplace à deux dimensions en coordonnées cartésiennes avec conditions aux limites sur un 
rectangle par la méthode de séparation des variables 

Exemples divers 

Rectangle soumis à des conditions de Dirichlet de valeurs fixées 

Rectangle soumis à des conditions de Dirichlet et de Neumann de valeurs fixées 

Rectangle soumis à des conditions de Dirichlet et de Robin de valeurs fixées 

Rectangle soumis à des conditions homogènes en y mixte de Robin/Dirichlet ou Robin/Neumann 

Rectangle soumis à des conditions homogènes mixte de Robin en y et inhomogène 

mixte en x 

Solution de l'équation de Laplace à deux dimensions sur des corps de dimension infini, 

Exemple trivial : un ruban infini entre x=0, x=Ix, x=0, condition aux limites homogènes de Dirichlet 


Solution de l'équation de Laplace à deux dimensions en coordonnées cylindrique (r, È) ou (r,z) avec conditions aux 
limites Dirichlet, Neumann ou Robin par la méthode de séparation des variables 
Singularité du système de coordonnées cylindriques 2D-3D et régularité du gradient 
équation de Laplace en coordonnées cylindrique polaire (r, 0) 
équation de Laplace en coordonnées cylindrique axiale (r,z) 
Cylindre plein (rz) soumis à des conditions aux limites de Dirichlet homogènes en z 
Cylindre creux (rz) soumis à des conditions aux limites de Dirichlet homogènes en z 
Cylindre plein (r,z) soumis à des conditions aux limites mixte de Dirichlet et Neumann homogènes en z 
Cylindre creux (rz) soumis à des conditions aux limites de Dirichlet homogènes en r 
Cylindre plein (rz) soumis à des conditions aux limites de Dirichlet homogènes en r 
Cylindre plein (rz) soumis à des conditions aux limites mixtes de Dirichlet et Neumann homogènes en r 
Cylindre plein (rz) soumis à des conditions aux limites mixtes de Dirichlet et Robin homogènes en r 
Cylindre creux (rz) soumis à des conditions aux limites mixtes de Robin homogènes en r 
Cylindre plein (rz) soumis à des conditions aux limites mixtes de Robin homogènes en z 
Cylindre creux (rz) soumis à des conditions aux limites mixtes de Robin homogènes en z 
Disque creux (r,0) soumis à des conditions aux limites mixtes de Dirichlet et Neumann homogènes en r 
Disque creux (r,0) soumis à des conditions aux limites de Dirichlet en r 
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Disque plein (r,0) soumis à des conditions aux limites de Dirichlet en r 
Disque plein (r,0) soumis à des conditions aux limites de Robin en r 
Disque creux (r,0) soumis à des conditions aux limites de Dirichlet en r 


Étude des problèmes aux limites de l'équation de Laplace sur des secteurs de disque 
Secteur de disque plein (r,0) soumis à des conditions aux limites de Dirichlet homogènes en 80 
Quelques cas simples de solutions par séparation des variables aux problèmes aux limites inhomogènes angulaires 
Secteur de disque plein (r), illimité dans le plan, soumis à des conditions aux limites de Dirichlet. 
Secteur de disque plein (r) soumis à des conditions aux limites mixtes de Robin homogènes en ÿ 
Secteur de disque annulaire (r,9) soumis à des conditions aux limites de Dirichlet homogènes en 8 
Secteur de disque plein (r,0) soumis à des conditions aux limites mixtes de Dirichlet, Neumann et Robin homogènes 
en 6 
Secteur de disque annulaire (r0) soumis à des conditions aux limites mixtes de Dirichlet, Neumann et Robin 
homogènes en 6 
Secteur de disque annulaire (r9) soumis à des conditions aux limites mixtes de Robin homogènes en ©, 
inhomogènes en r 
Secteur de disque annulaire (r) soumis à des conditions aux limites mixtes de Robin inhomogènes en r, homogènes 
en Ÿ 
Problèmes homogènes dans la dimension radiale pour un secteur de disque annulaire 
Généralisation des conditions aux limites homogènes radiales 
cas général secteur annulaire avec conditions aux limites homogènes mixtes de Robin 
Extension au cas cartésien avec l'angle d'ouverture du cône a x 


Étude des problèmes aux limites de l'équation de Laplace en configuration cartésiennes 3D 


Étude des problèmes aux limites de l'équation de Laplace en configuration cylindrique polaire et axiale trois 
dimensions, cas du cylindre complet plein ou creux, de hauteur z, cas d'un secteur de cylindre plein ou creux de 
hauteur z 
Singularité du système de coordonnées cylindriques 3D et régularité du gradient 
Orthogonalité bidimensionnelle des fonctions propres sur un problème homogène cylindrique dans les deux 
dimensions radiale r et angulaire Ÿ 
Le cylindre plein complet 
Problème homogène dans la dimension axiale z 
Problème homogène dans la dimension radiale r 
Le cylindre creux complet 
Problème homogène axial sur le cylindre complet creux 
Problème homogène radial sur le cylindre complet creux 
Secteur cylindrique plein 
Solution du problème inhomogène axial pour le secteur de cylindre plein: 
Solution du problème inhomogène radial pour le secteur de cylindre plein 
Solution du problème inhomogène angulaire pour le secteur de cylindre plein 
Secteur cylindrique creux 
Solution du problème inhomogène axial pour le secteur de cylindre creux 
Solution du problème inhomogène radial pour le secteur de cylindre creux 
Solution du problème inhomogène angulaire pour le secteur de cylindre creux 
Fonctions de Bessel modifiées d'ordre purement imaginaire 
Fonctions de Bessel d'ordre purement imaginaire 
Quelques intégrales des fonctions de Bessel 
Les formules les plus couramment utilisées des intégrales des fonctions de Bessel 
Utilisation des fonctions hypergéométriques 
Utilisation des fonctions de Lommel (cas particuliers des fonctions de Struve) 
Calcul des intégrales des fonctions de Bessel à l'aide des fonctions de Lommel 
Utilisation des fonctions de Struve 


Méthode de résolution des problèmes aux limites de l'équation de Laplace en coordonnées cylindriques 3D sur des 
domaines non bornés par les transformations intégrales de Hankel 
Théorème de Hankel avec les fonctions de Bessel de première espèce 
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Transformation de Hankel avec les fonctions de Bessel de deuxième espèce 

Transformation sur l'intervalle [x0,+c°) appelée transformation de Weber 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème de Dirichlet sur le demi-espace z>0 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème mixte de Robin sur le demi-espace z>0 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème de Dirichlet sur une bande centrée de hauteur 2z0 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème de Dirichlet électrostatique dans l'espace plan de hauteur 220, potentiel dû à une charge placée à 
l'origine des coordonnées 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème de Neumann électrostatique dans de demi espace z>0, potentiel dû à un disque de rayon rO posé à la 
terre appliquant un courant continu 


Solution de l'équation de Laplace en coordonnées sphériques 2D, 3D avec conditions aux limites Dirichlet, Neumann 
ou Robin par la méthode de séparation des variables 
Singularité du système de coordonnées sphérique 3D et régularité du gradient 
Problème intérieur et extérieur, sphère pleine (r,0,@) soumise à des conditions aux limites de Dirichlet inhomogènes 
en r=lr. ne dépendant que de l'angle 8 et pas de l'angle azimutal @ 
Application à un problème d'électrostatique 
Fonction génératrice des polynômes de Legendre et des fonctions associées de Legendre 
Problème intérieur, Sphère pleine (r,90,) soumise à des conditions aux limites de Dirichlet inhomogènes en r=lr. ne 
dépendant que de l'angle @ 
Sphère pleine (r80) soumise à des conditions aux limites de Neumann inhomogènes en r=lr 
Sphère pleine (r0) soumise à des conditions aux limites de Robin inhomogènes en r=lr. 
Hémisphère pleine (0) soumise à des conditions aux limites de Dirichlet inhomogènes 
Sphère creuse (r,0) soumise à des conditions aux limites de Dirichlet inhomogènes en r=lr2 et Dirichlet homogènes 
en r=lr1. 
Sphère creuse (r,9) soumise à des conditions aux limites de Robin inhomogènes en r=lr2 et Dirichlet homogènes en 
r=lr1. 
Hémisphère creux (r0) soumise à des conditions aux limites de Dirichlet inhomogènes en r=lr2 et homogènes en 
r=lr1. 
Cône sphérique plein d'angle 00, en coordonnées sphériques (r0) soumis à des conditions aux limites de Dirichlet 
inhomogènes en r=lr. homogènes en 00 
Cône sphérique plein d'angle 00, en coordonnées sphériques (r90) soumis à des conditions aux limites de Dirichlet 
inhomogènes en r=lr. et de Neumann homogènes en 600 
Considération sur les valeurs propres du problème aux limites solution d'équation transcendantales. 
Estimation des valeurs propres du problème de Dirichlet sur le cône sphérique 
Valeurs des zéros des fonctions associées de Legendre pour m> 0 
Estimation des racines pour quelques valeurs limites de 90 : 90=x/2 
Estimation des racines pour quelques valeurs limites de 90 : 90 proche de zéro 
Estimation des racines pour quelques valeurs limites de 90 : 90 proche de n 
Valeurs propres du problème de Neumann sur le cône sphérique pour des valeurs moyennes de l'angle d'ouverture 
n/6<00<5n/6 
Problème aux limites mixte de Robin en configuration conique-sphérique 
Problèmes aux limites sur des sections coniques-sphériques 
Quelques propriétés des fonctions de Legendre de première et deuxième espèce pour les degrés Àn non entiers 
pour des problèmes aux limites de Dirichlet 
Quelques propriétés des fonctions de Legendre de première et deuxième espèce pour les degrés Àn non entiers 
pour des problèmes aux limites de Neumann 
Quelques propriétés des fonctions de Legendre de première et deuxième espèce pour les degrés Àn non entiers 
pour des problèmes aux limites mixtes > de Robin dans le cas d'un problème entièrement symétrique 
Calcul de la norme des fonctions propres du problème aux limites de Dirichlet 
Calcul de la norme des fonctions propres du problème aux limites de Neumann 
Calcul de la norme des fonctions propres du problème aux limites de Robin 
Extension aux cas de l'hémisphère porté a la condition de Dirichlet homogène nulle sur la base 
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Tableau synoptique des normes des fonctions propres angulaires d'un problème aux limites sur une section 
sphériques 
Solutions de l'équation de Laplace à l'extérieur d'un domaine Q en coordonnées sphériques 
Solution des problèmes extérieur et intérieur électrostatique sur des corps géométrique de révolution connaissant la 
valeur du potentiel sur l'axe de révolution (d'après Jackson, Classical Electrodynamics) 
Le problème des deux sphères de rayon de potentiel inégal placées à distance l'une de l'autre ou au dessus d'un 
plan conducteur 
Problèmes aux limites de Laplace homogène en dimension radiale sur des sections coniques pleines et creuses 
Orthogonalité et normalisation des fonctions propres radiales construites 
Calcul des normes en théorie de Sturm-Liouville 
Formes des fonctions angulaires associées 
Problème aux limites sur une section conique-sphérique creuse avec des conditions aux limites en dépendance 
azimutale 
Fonctions propres radiales 
Les fonctions angulaires azimutales 
Les fonctions angulaires en ÿ : 
Représentation intégrale de la solution du problème de Dirichlet inhomogène sur la surface d'une section conique- 
sphérique a trois dimensions 
intégrales et représentations intégrales des fonctions coniques de Mehler d'ordre O 
Autres représentations intégrales des fonctions coniques de Mehler d'ordre supérieur 
Exemples divers de représentations intégrales 


Sphère et hypersphère 
Problèmes aux limites unidimensionnels en dimension radiale 
Solutions de l'équation de Laplace à l'extérieur d'un domaine Q ultra-sphérique 
Problèmes aux limites de Laplace sur une hypersphère avec dépendance à une des co-latitudes 
Solution des problèmes aux limites homogènes en dimension radiale 
Solution des problèmes aux limites homogènes en dimension angulaire ou << inhomogènes en dimension 
radiale sur l'intervalle angulaire complet [0,n] > 
Construction des fonctions de Gegenbauer de première et deuxième espèce sur l'intervalle [-1,1] 
Construction d'une seconde solution indépendante de l'équation de Gegenbauer sur l'intervalle1<z< co 
Wronskien sur l'intervalle ]-1,1[ 
Wronskien pour z > 1 
Lien des polynômes de Gegenbauer avec les polynômes de Jacobi 
Orthogonalite des fonctions de Gegenbauer de première et de deuxième espèce 
Formule de Rodrigues et fonctions génératrices des polynômes de Gegenbauer 
Formules de récurrence des fonctions de première espèce et des polynômes de Gegenbauer 
Formules de dérivation des fonctions de Gegenbauer de première espèce 
Lien avec d'autres polynômes dans la littérature scientifique : 
Calcul de quelques intégrales indéfinies des fonctions de Gegenbauer de première et deuxième espèce 
Intégrales indéfinies des fonctions de Gegenbauer de première et deuxième espèce pondérées 
intégrales indéfinies des fonctions de Gegenbauer a l'aide d'équations conjuguées 
intégrales indéfinies des fonctions de Gegenbauer a l'aide de fonctions h(x) quelconques 
Calcul des normes des fonctions de Gegenbauer, des fonctions propres d'un problème aux limites, approche de 
Sturm-Liouville 
cône de révolution ultra-sphérique 
Problème aux limites a l'intérieur ou l'extérieur d'une hyper-sphère 
problème aux limites intérieur sur une hyper-sphère soumise à des conditions aux limites inhomogènes de 
Dirichlet dépendante uniquement de l'angle ĝn 
problème aux limites extérieur sur une hyper-sphère soumise à des conditions aux limites inhomogènes de 
Dirichlet dépendante uniquement de l'angle n 
Application a un problème d'électrostatique a N-dimensions 
problème aux limites extérieur et intérieur sur une hyper-sphère soumise à des conditions aux limites 
inhomogènes de Robin dépendante uniquement de l'angle ÿn 
problème aux limites intérieur sur une hyper-sphère creuse soumise à des conditions aux limites 
inhomogènes de Robin dépendante uniquement de l'angle n. de part et d'autre. 
Construction des fonctions de Gegenbauer de première et deuxième espèce, formules de liaison 
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Construction des solutions des problèmes aux limites et parité des fonctions propres 
Cas des fonctions propres solutions du problème aux limites homogène de Dirichlet sur une section conique 
hyper-sphérique 
Cas des fonctions propres solutions du problème aux limites homogène de Neumann sur une section conique 
hyper-sphérique 
problème aux limites a l'intérieur d'un cône hyper-sphérique d'angle ouvert 90 fini 
problème aux limites a l'intérieur d'un cône hyper-sphérique d'angle ouvert 00 non bornée radialement 
Valeur des dérivées paramétriques des fonctions de Gegenbauer dans le cas d'un cône d'angle d'ouverture 
Ÿ0=71/2 
problème aux limites a l'intérieur d'un cône hyper-sphérique d'angle ouvert 90 bornée ou non radialement, 
condition aux limites constantes. 
problème aux limites a l'intérieur d'un cône hyper-sphérique d'angle ouvert 90 non borne radialement avec des 
conditions inhomogènes de Dirichlet sur la surface latérale du cône 


Exemple de valeurs propres d'un problème aux limites sur un cône ultra-sphérique d'angle ouvert 90 ou une section 
conique sphérique d'angle 1 et ÿ2 
Extension aux cas de l'hémisphère a N-dimensions porte a la condition de Dirichlet homogène nulle sur la base 
Solution des problèmes aux limites homogènes en dimension radiale sur un cône ultra-sphérique 
Calcul des normes en théorie de Sturm-Liouville 
Formes des fonctions angulaires associées 
Hypothèse sur les formes des fonctions angulaires de valeurs propres nulle suivant les valeurs du paramètre 
dimensionnel À 
Remarque sur d'autre choix de seconde solution 


Représentation intégrale de la solution du problème de Dirichlet inhomogène sur la surface d'un cône ultra- 
sphérique 
Comportement des solutions à la singularité 9 =0 


Représentation des potentiels sous forme intégrale 


Quelques considérations sur les dérivées des fonctions de Gegenbauer et Legendre de première espèce par rapport 
à son degré. Valeurs de la dérivée aux degrés entiers pour des ordres entiers et demi-entiers 

Dérivées premières par rapport au degré, fonctions de Legendre de degré quelconque 

Dérivées premières par rapport au degré, fonctions de Legendre de degré entier 


Dérivées paramétriques des polynômes de Gegenbauer d'ordre demi-entier (à finir mettre dans le texte) 
Expression formelle des dérivées paramétriques selon le degré pour les fonctions/polynômes de Gegenbauer, 
valeur au degré 0 
Expression formelle des dérivées paramétriques selon le degré pour les fonctions/polynômes de Gegenbauer, 
valeur au degré 1 
Formule de récurrence sur les dérivées paramétriques des polynômes de Legendre, demonstration par 
récurrence de la formule de Bromwhich pour les polynômes de Legendre 
Formule de récurrence sur les dérivées paramétriques des polynômes de Gegenbauer de degré demi-entier, 
nouvelle formule de Bromwhich pour les polynômes de Gegenbauer de degré demi-entier 
Application de la formule de Bromwhich au calcul des fonctions de Gegenbauer de deuxième espèce de 
degré entier et d'ordre demi-entier 
Table des nouveaux polynômes ultra-sphériques de Bromwhich pour des ordres demi-entiers 
Applications aux normes de fonctions propres de problèmes aux limites en physique mathématique 


Dérivées paramétriques des polynômes de Gegenbauer d'ordre entier 
Application au calcul des fonctions de Gegenbauer de deuxième espèce de degré entier et d'ordre entier 
Récurrence des dérivées paramétriques des polynômes de Gegenbauer sur les ordres entiers ou demi-entiers 
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Hypothèse de récurrence des dérivées paramétriques des polynômes de Gegenbauer pour les ordres entiers 
Application au calcul des fonctions de Gegenbauer de deuxième espèce degrés quelconque et d'ordre entier 
Table des polynômes ultra-sphériques pour des ordres entiers 

Application aux calculs des normes des fonctions propres de problème aux limites ultra-sphériques 


Problèmes aux limites de Laplace sur des systèmes de coordonnées curvilignes orthogonales, engendrant des 
surfaces de révolution 
Séparabilité simple de l'équation de Helmholtz ou de l'équation de Laplace 
R-Séparabilité de l'équation de Laplace 
équation différentielle des fonctions de révolution 
Différentes solutions de l'équation différentielle pour différents systèmes de coordonnées 
Application a des systèmes de révolution connus, simple séparabilité 
Le système cylindrique 
Le système sphérique 
Le système sphéroïdal allongé 
Le système sphéroïidal aplati 
Le système parabolique de révolution ou paraboloïdal 


Application à des systèmes de révolution connus, R séparabilité 
Les coordonnées toroidales 
Les coordonnées de type toroïdal, contribution de René Lagrange < Les familles de surface de révolution qui 
possèdent des harmoniques, 1939 » 
Les coordonnées bi-sphériques 
Les coordonnées inverses sphéroïdales allongées 
Les coordonnées inverses sphéroïdales aplaties 
Le système cardioïde de révolution ou parabolique de révolution inverse > 
Le système sphère-tangentielle de révolution ou < cylindrique inverse > 
Les coordonnées cyclidiques, contribution de René Lagrange < Les familles de surface de révolution qui 
possèdent des harmoniques, 1939> 
Quelques formules utiles 


Solution de l'équation de Laplace en coordonnées sphéroïdales allongées (prolate sphéroïdal) ou aplaties (oblate 
sphéroïdal) à trois dimensions avec conditions aux limites Dirichlet, Neumann ou Robin par la méthode de 
séparation des variables, passage à deux dimensions 


Le système de coordonnées sphéroïdales allongées (prolate sphéroïdal) 


Iso-surfaces du système de coordonnées 
Singularité du système de coordonnées sphéroïdales allongées et régularité du gradient 
Solutions de l'équation de Laplace par séparation des variables 
Exemples divers sur des sphéroïdes et “"Hémisphéroïdes" allongés pleins ou creux 
Exemples divers sur des cônes sphéroïdaux ou sections coniques sphéroïdales allongés plein ou creux (domaines 
bordés par des hyperboloides à deux nappes, tronqués), soumis à des conditions aux limites homogènes 
angulaires et inhomogènes radiales 
Exemple : problème de Dirichlet sur un cône sphéroïdal non borné ou hyperboloïde à deux nappes infini 
Problèmes aux limites sur des hyperboloïides à deux nappes (cônes sphéroïdaux ou section coniques- 
sphéroidales) homogènes dans la dimension radiale n, sans dépendance azimutale 
Représentation intégrale des solutions d'un problème aux limites inhomogène angulaire sur l'hyperboloïde de 
révolution à deux nappes infini 

Théorème de Mehler-Fock et transformation intégrale de Mehler-Fock 

Transformation intégrale de Mehler-Fock d'ordre supérieur 

Solution du problème aux limites de Dirichlet sur l'intervalle [0,+c°[ 
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Solution du problème aux limites de Dirichlet sur l'intervalle [no,+c°[ 

Représentation intégrale modifiée dans le cas d'un problème aux limites sur [n0,+ce[, première forme 
d'après l'article de N.A.BELOVA.la.S.UFLIAND, « Dirichlet problem for a toroidal segment », PMM Vol. 31, 
No. 1, 1967, pp. 59-63 » 

Représentation intégrale modifiée dans le cas d'un problème aux limites sur [n0,+ce[, première forme 
d'après l'article de A.F.Ulitko « On a generalization of the integral transformation of Mehler-Fock », SOVIET 
APPLIED MECHANICS, Prikladnaya Mekhmnika, Vot. 3, No. 5, pp. 45-49, 1967» 

Équivalence des deux représentations intégrales de Belova-Ufliand et d'Ultiko 

Forme intégrale des solutions des problèmes aux limites sur des hyperboloïdes tronqués 


Le système de coordonnées sphéroïdales aplaties (oblate sphéroïdal) 


Iso-surfaces du système de coordonnées 
Singularité du système de coordonnées sphéroïdales aplaties et régularité du gradient 
Solutions de l'équation de Laplace par séparation des variables 
Exemples divers sur des sphéroïdes et “"hémisphéroïdes" aplatis pleins ou creux 
Exemples divers sur des cônes sphéroïdaux ou sections coniques sphéroïdales allongés plein ou creux (domaines 
bordés par des hyperboloïides à une nappe, tronqués), soumis à des conditions aux limites homogènes 
angulaires et inhomogènes radiales 
Tableau synoptique des normes des fonctions propres angulaires d'un problème aux limites sur une section 
sphérique 
Problèmes aux limites sur des hyperboloïdes à une nappe, homogènes dans la dimension radiale n, sans 
dépendance azimutale 
Considérations générales 
Quelques propriétés des fonctions coniques et fonctions coniques associées de Mehler pour un argument 
purement imaginaire 
Représentation intégrale des solutions d'un problème aux limites inhomogène angulaire sur l'hyperboloïde 
de révolution à une nappe 
Transformation intégrale modifiée de Mehler-Fock 
Exemples de développement intégrale d'une fonction limite par des fonctions coniques de Mehler 
d'argument imaginaire 
Applications de la transformée intégrale modifiée de Mehler-Fock à des problèmes aux limites sur des 
hyperboloïdes à une nappe 


Solution de l'équation de Laplace à deux dimensions en coordonnées elliptiques cylindrique et polaire avec 
conditions aux limites Dirichlet, Neumann ou Robin par la méthode de séparation des variables 


Le système elliptique polaire ou cylindrique (2D, 3D) 

Relations liant une ellipse de grand axe a et petit axe b, d'excentricité e, avec le système de coordonnées elliptiques 
Jacobien du système de coordonnées 

Singularité du système de coordonnées elliptiques cylindriques 3D et régularité du gradient 

Singularité du système de coordonnées elliptiques cylindriques 2D et régularité du gradient 

Conditions de continuité et dérivabilité sur la ligne interfocale pour des solutions admissibles à 2 et 3 dimensions 
Séparation des variables en coordonnées elliptiques polaires 2D 

Séparation des variables en coordonnées elliptiques cylindriques 3D 

Problème aux valeurs propres et formes des solutions 

Périodicité des solutions et conséquence de la symétrie des conditions aux limites sur la forme des solutions 
Formulaires des fonctions de Mathieu et fonctions modifiées de Mathieu de première et deuxième espèces 
Méthode de calcul des coefficients des fonctions de Mathieu angulaires 

Fonctions de Mathieu de première espèce (partie angulaire), cas q>0 : 

Propriétés de parité, de périodicité, de translation des fonctions de Mathieu de première espèce (partie angulaire), 
cas q>0 

Fonctions de Mathieu de première espèce (partie angulaire), cas q<0 : 

Fonctions de Mathieu de deuxième espèce (partie angulaire), cas q<0 : 

Changement de signe dans les coefficients du développement de Fourier 

Fonctions modifiées de Mathieu de première espèce (partie radiale), cas q>0 
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Fonctions modifiées de Mathieu de deuxième espèce (partie radiale), cas q>0 

Fonctions modifiées de Mathieu de première espèce (partie radiale), cas q<0 

Fonctions modifiées de Mathieu de deuxième espèce (partie radiale), cas q<0 

Dérivées des fonctions modifiées de Mathieu de première espèce (partie radiale), cas q>0 

Dérivées des fonctions modifiées de Mathieu de deuxième espèce (partie radiale), cas q>0 

Dérivées des fonctions modifiées de Mathieu de deuxième espèce (partie radiale), cas q<0 

Dérivées des fonctions modifiées de Mathieu de première espèce (partie radiale), cas q<0 

Relation d'orthogonalité unidimensionnelle des fonctions de Mathieu 

Relation d'orthogonalité bidimensionnelle des fonctions de Mathieu 

Évaluation de certaines intégrales des fonctions de Mathieu à partir de leur décomposition 

Dégénérescence de l'ellipse en cercle et formes des solutions radiales et angulaires à trois dimensions 
Dégénérescence de l'ellipse en cercle et formes des solutions radiales et angulaires à deux dimensions 

Exemples de solutions de l'équation de Laplace en coordonnées elliptiques cylindriques à deux dimensions 

Exemples à deux dimensions : problème intérieur de Dirichlet sur une section elliptique pleine et bornée avec une 
fonction limite quelconque 

Exemple à deux dimensions : problème intérieur de Dirichlet sur une section elliptique bornée creuse avec une 
fonction limite quelconque, condition aux limites homogène sur le rayon intérieur 

Exemples à deux dimensions : problème de Dirichlet sur une bande elliptique infini, comprise entre deux hyperboles 
Exemple à trois dimensions : problème de Dirichlet sur une section elliptique pleine et bornée avec une fonction 
limite quelconque, et une hauteur tranchée en z avec des conditions homogènes sur les sections haute et basse 
Exemple à trois dimensions : problème de Dirichlet sur une section elliptique creuse et bornée avec une fonction 
limite quelconque, et une hauteur tranchée en z avec des conditions homogènes sur les sections haute et basse et 
dans la paroi latérale intérieure 

Exemple à trois dimensions : problème aux limites sur une section elliptique pleine et bornée avec une condition aux 
limites homogène de Dirichlet sur la face cylindrique, et une hauteur tranchée en z avec des conditions 
inhomogènes sur la section haute et homogène sur la section basse (à changer dans le texte) 

Exemple à trois dimensions : problème aux limites sur une section elliptique pleine et bornée avec une condition aux 
limites homogène de Neumann sur la face cylindrique, et une hauteur tranchée en z avec des conditions 
inhomogènes sur la section haute et homogène sur la section basse.(à changer dans le texte) 


Solution de l'équation de Laplace à deux dimensions en coordonnées paraboliques polaires 2D (u,v) et cylindriques 
3D (u,v,z) avec conditions aux limites Dirichlet, Neumann ou Robin par la méthode de séparation des variables 


Le système parabolique cylindrique (u,v,z) (2D,3D) 

Propriétés des fonctions de Weber cylindriques paires We(p,z), impaires Wo(p,z), fonctions U(a,z), V(a,z), W{a,z) et 

Dp(z) 

Propriétés des fonctions paraboliques cylindriques de Weber 

Propriétés des fonctions paraboliques cylindriques U(a,z) et V(a,z) 

Propriétés des fonctions paraboliques cylindriques Dp(z) de Whittaker-Watson 

Lien avec les fonctions de Hermite-Gauss, fonction génératrice, intégrale de Fourier 

Singularité du système de coordonnées paraboliques cylindrique 3D et régularité du gradient 

Régularité du gradient avec les fonctions paraboliques cylindriques de Weber 

Régularité du gradient avec les fonctions paraboliques cylindriques D 

Régularité du gradient avec les fonctions paraboliques cylindriques U(a,z) et V(a,z) 

Régularité du gradient avec les fonctions de Bessel solutions des cas n°5 et 6 

Singularité du système de coordonnées parabolique cylindrique en 2D, régularité du gradient à 

deux dimensions 

Orthogonalité des solutions de l'équation de Laplace en coordonnées paraboliques cylindriques 

Problème aux limites à deux dimensions en coordonnées paraboliques cylindriques 
Exemple à deux dimensions : problème de Dirichlet sur une section parabolique pleine et bornée avec une 
fonction limite quelconque. 
Exemple à deux dimensions : problème d'une bande infinie soumise à des conditions homogène de Dirichlet en 
v1 et constante en v2 
Exemple : soit en coordonnées parabolique 2D le problème d'un demi-espace parabolique soumis à des 
conditions inhomogènes de Dirichlet quelconque en vO, solution intégrale. 
Exemple : soit en coordonnées parabolique 2D le problème d'un demi-espace parabolique soumis à des 
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conditions inhomogènes de Dirichlet quelconque en 1=u0 et 1=-u0, solution 
intégrale. 


Problème aux limites à trois dimensions en coordonnées paraboliques cylindriques 
Exemple : soit en coordonnées parabolique 3D le problème de Dirichlet sur un domaine parabolique cylindrique 
de dimension finie et de configuration géométrique symétrique dans les coordonnées proprement parabolique, 
et homogène à la surface latérale du cylindre 
Exemple : problème intérieur convexe inhomogène de Dirichlet en v=vO , dans un domaine parabolique non 
borné en L, tranché en z de hauteur z finie , avec condition homogène de Dirichlet en z. 
Exemple : problème intérieur convexe inhomogène de Dirichlet en v=vO dans un domaine parabolique non borné 
en u, tranché en z de hauteur z finie, avec condition homogène de Neumann en z. 
Exemple : Problème extérieur concave, inhomogène de Dirichlet en v=v0 dans un domaine parabolique non 
borné en u, tranché en z de hauteur z finie, avec condition homogène de Dirichlet en z 
Exemple : Problème extérieur concave inhomogène de Dirichlet en v=vO dans un domaine parabolique non 
borné en L, tranché en z de hauteur z finie, avec condition homogène de Neumann en z 
Exemple : problème intérieur convexe, inhomogène de Dirichlet en v=v0 dans un domaine parabolique non 
borné en L, de hauteur infinie en z 
Exemple : problème extérieur concave, inhomogène de Dirichlet en v=vO , dans un domaine parabolique, de 
hauteur z infinie. 
Exemple : Problème intérieur de Dirichlet dans un cylindre de sections paraboliques, de hauteur z finie, borné en 
u et v, -u0 < u <u0 ,0<v <v0 , homogène en u, inhomogène en v 
Exemple de valeurs propres dans une configuration de section parabolique symétrique de largeur Ix pour le 
problème aux limites paires 
Exemple de valeurs propres dans une configuration de section parabolique symétrique de largeur Ix pour le 
problème aux limites impaires 
Exemple : Problème intérieur de Dirichlet dans un cylindre de sections paraboliques, de hauteur z finie, borné en 
u et v, -u0 < u <u0 ,0<v <v0 , homogène en v, inhomogène en u 
Exemple de valeurs propres pour le problème aux limites paires 
Exemple de valeurs propres pour le problème aux limites impaires 


Solution de l'équation de Laplace en coordonnées toroïdales (n,9,@) avec conditions aux limites Dirichlet, Neumann 
ou Robin par la méthode de séparation des variables 


Le système toroïdal (n,9,¢) (3D) 
Iso-surface du système de coordonnées toroïdales 
Valeurs particulières du système toroïdal (à changer dans le texte) 
Comportement des iso-surfaces Ÿ=00 
Comportement des iso-surfaces n=n0 
Singularité du système de coordonnées toroïdales et régularité du gradient 
Forme traditionnelle de la R-Séparabilité en coordonnées toroïdales (R-Séparabilité-1) 
Nouvelle forme de la R-Séparabilité en coordonnées toroïdales (R-Séparabilité-2) 
Contraintes et conditions imposées pour les solutions de l'équation de Laplace pour les problèmes aux limites dans 
le système de coordonnées toroïdales 
Problème intérieur de Dirichlet à deux et trois dimensions 
R-séparabilité 1 : problème intérieur dépendant des variables n et Ÿ 
Problème intérieur de Dirichlet sur un tore avec une condition aux limites impaire en ÿ 
Solution général au problème de Dirichlet sur un-tore complet pour des conditions aux limites en admettant 
un développement en série de Fourier 
Extension au problème d'un demi-tore dont la base est portée à la valeur 0 
Problème intérieur de Dirichlet sur un tore complet et creux 
Problème intérieur homogène et inhomogène de Dirichlet sur un tore complet et creux 
Problème intérieur homogène et inhomogène de Dirichlet sur un demi tore complet et creux 
Problème intérieur homogène de Neumann et inhomogène de Dirichlet sur un tore complet et creux 
Divers exemples tirés de la littérature 
Homogénéisation de l'équation de poisson pour un système de coordonnées orthogonales de révolution 
Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique 
uniforme le long de l'axe z, porté au potentiel 0, exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, 
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Y.S.Ufliand « Problems of Mathematical physics » (problème n° 494 page 236) 

Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique uniforme le long 
de l'axe z, porté au potentiel VO 

Exemple d'électrostatique : problème extérieur de Dirichlet, tore dans un champ électrique uniforme le long 
de l'axe x ou y, porté au potentiel O 

Exemple d'électrostatique : problème extérieur de Dirichlet, potentiel extérieur à la surface d'un tore 
conducteur parfait induit par une charge placée à l'origine des coordonnées et donc porté au potentiel 0. 
Exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, Y.S.Ufliand « Problems of Mathematical physics » 
(problème n° 495 page 237) 

Exemple régime stationnaire équation de la chaleur en présence de source de chaleur uniforme: problème 
intérieur de Dirichlet, à l'intérieur d'un tore dont la surface est à une valeur nulle avec une répartition 
uniforme de source de chaleur. Exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, Y.S.Ufliand « 
Problems of Mathematical physics » (problème n° 498 page 238) 

Exemple gravitationnel: problème extérieur de Dirichlet, champ gravitationnel d'un tore de densité massique 
constante p. Exemple tiré de l'ouvrage de N.N.Lebedev, I.P.Skalkaya, Y.S.Ufliand « Problems of Mathematical 
physics » (problème n° 499 page 238) 


Solution sur une lentille sphérique symétrique assimilable à un problème homogène angulaire Problème 
intérieur de Dirichlet à trois dimensions avec dépendance azimutale : équivalence des solutions en R- 
Séparabilité-1 et R-Séparabilité-2 
R-séparabilité 2 : problème intérieur dépendant des variables n, Ÿ et D 
Problème extérieur de Dirichlet 

R-séparabilité 1 problème extérieur de Dirichlet 

R-séparabilité 2 : problème extérieur dépendant des variables n et D 
Remarques sur les fonctions propres sur des tores complets pleins en fonction de la décomposition par symétrie 
des conditions aux limites 
Remarque sur la solution du problème aux limites inhomogène de Neumann ou mixte de Robin Solutions des 
problèmes aux limites homogènes angulaires Ÿ sur des sections « coniques » 
toroïdales 

Quelques illustrations géométriques des domaines à l'étude 

Solutions sur les domaines radiaux n €[In0 ,+ ce ] 

Solutions sur les domaines radiaux n €[0,In0] 

Solutions sur les domaines radiaux n €fin1,1n2] 


Solutions des problèmes aux limites inhomogènes angulaires Ê sur des sections « coniques » toroïdales, 
conditions homogène radiales 
Discussion sur le cas de l'intervalle n€[0,1n0], z=Cosh(n)->z£[1,z0] (voir système sphéroïdal) 
Discussion sur le cas de l'intervalle n€[ In1 , In2 ], z=Cosh(n) -> z£€[z1=Cosh( In1 ), z2=Cosh( In2 ) ] (voir 
système sphéroïdal) 
Discussion sur le cas de l'intervalle n€[ln0 ,+ ce ], z=Cosh(n)->z€[z0=Cosh(In0 ),+ ce ] (voir système sphéroïdal) 
Discussion sur le cas de l'intervalle n€[0,+ ce ] (voir système sphéroïdal) 
Problème aux limites sur l'intervalle [0,1In0] 
Problème aux limites sur l'intervalle [In0 ,+ ce ] 
Problème aux limites sur l'intervalle [In1,1n2] 


Solutions des problèmes aux limites sur des domaines union, intersection et différence de sphères, lentilles 
sphériques 
Modélisation d'intersections de sphères (lentilles sphérique) 
Modélisation d'unions de sphères 
Modélisation de différences de sphères 
Exemple : solution intégrale du problème de Dirichlet sur le domaine intérieur de la différence de deux sphères 
par la transformation de Mehler-Fock (à changer dans le texte) 
Exemple : solution du problème de Dirichlet sur le domaine intérieur d'un hémisphère plein 
Exemple : solution du problème de Dirichlet sur le domaine extérieur de la différence de deux sphères 
Exemple : solution du problème de Dirichlet sur le domaine extérieur d'un hémisphère plein 
Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur de la différence de sphères : 
potentiel dû à une charge placée à l'origine des coordonnées 
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Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur de la différence de sphères : 
potentiel extérieur d'un corps au potentiel nul et placé dans un champ électrostatique uniforme orienté selon 
l'axe z négatif 

Exemple : solution du problème de Dirichlet sur le domaine extérieur d'une lentille sphérique, ou également 
d'une union sphérique 

Exemple : solution du problème de Dirichlet sur le domaine intérieur d'une lentille ou d'une union sphérique 
Exemple : solution du problème de Dirichlet sur le domaine extérieur d'une coupole ou coquille sphérique 
Formulaire de passage entre les coordonnées toroïdales et les coordonnées cylindriques 3D 

Exemple : solution du problème de Dirichlet électrostatique sur le domaine extérieur d'une coupole sphérique : 
potentiel dû à une charge placé à l'origine des coordonnées 


Représentations intégrales des fonctions de Legendre et des fonctions associées de Legendre Autres 
représentations intégrales des fonctions de Legendre de degré demi-entier Quelques développement en série 
utilisant les fonctions toroïdales, dans le cadre de l'inverse des distances en coordonnées cylindriques 3D 


Solution de l'équation de Laplace en coordonnées parabolique de révolution (u,v,@) avec conditions aux limites 
Dirichlet, Neumann ou Robin par la méthode de séparation des variables 


Le système parabolique de révolution (u,v,) (3D) 
Singularité du système de coordonnées parabolique de révolution 2D et 3D et régularité du gradient 
Séparation des variables dans le système parabolique de révolution (à changer dans le texte) 
Paraboloïde plein de révolution soumis à des conditions aux limites indépendantes de l'angle azimutal 
Exemples divers (à changer dans le texte) 
Exemple : problème aux limites de Dirichlet sur un paraboloïde infini 
Exemple : un autre problème aux limites de Dirichlet sur un paraboloïde infini 


Paraboloïde plein de révolution soumis à des conditions aux limites dépendantes de l'angle azimutal 
Exemples divers (à changer dans le texte) 


Paraboloïde creux de révolution soumis à des conditions aux limites dépendantes ou 
indépendantes de l'angle azimutal 
Problème homogène dans la direction u 
Problème homogène dans la direction v 
Prise en compte de la valeur propre nulle (Problème homogène dans la direction u et v) 
Problème de Dirichlet homogène dans la direction u avec symétrie azimutale Y+,X+ 
Problème de Dirichlet homogène dans la direction v avec symétrie azimutale Y+,X+ 
Problème de Dirichlet avec valeur en limite indépendante de l'angle azimutal 
Problème de Dirichlet avec valeur en limite constante 
Exemple : Problème homogène en u de Neumann, homogène de Dirichlet en v=v1, et inhomogène de Neumann 
v2 à flux constant 
Exemple : Problème homogène en u de Neumann, homogène de Dirichlet en v=v1, et inhomogène de Dirichlet 
v2 à valeur constante 


Secteur de paraboloïde plein de révolution soumis à des conditions aux limites dépendantes ou indépendantes de 
l'angle azimutal 

Le problème inhomogène en u sur le secteur de paraboloïde plein de révolution 

Le problème inhomogène en v sur le secteur de paraboloïde plein de révolution 

Le problème inhomogène en ® sur le secteur de paraboloïde plein de révolution 


Secteur de paraboloïde creux de révolution soumis à des conditions aux limites dépendantes ou indépendantes de 
l'angle azimutal 

Le problème inhomogène en u sur le secteur de paraboloïde creux de révolution 

Le problème inhomogène en v sur le secteur de paraboloïde creux de révolution 

Le problème inhomogène en ® sur le secteur de paraboloïde creux de révolution 


Domaine non bornée ayant comme frontière un paraboloïde de révolution soumis à des conditions aux limites 
inhomogènes dépendantes ou indépendantes de l'angle azimutal 


Problèmes aux limites de Laplace en 2D-3D - Considérations générales — p 12 


Théorème de Hankel avec les fonctions de Bessel de première espèce 

Transformation de Hankel avec les fonctions de Bessel de deuxième espèce 

Transformation sur l'intervalle [x0,+ + ) appelée transformation de Weber 

Application aux problèmes aux limites intérieur et extérieur de Dirichlet sans dépendance azimutal domaine 
frontière d'un paraboloïde complet u£€[0, ce] 

Application des transformées de Hankel à la résolution des problèmes aux limites de l'équation de Laplace : 
problème de Dirichlet électrostatique dans l'espace intérieur d'un paraboloïde v=v 0, potentiel dû à une charge 
placée au foyer de la parabole 

Application aux problèmes aux limites inhomogène en v intérieur et extérieur de Dirichlet avec dépendance 
azimutal 

Application aux problèmes aux limites intérieur et extérieur inhomogène de Neumann en v sans dépendance 
azimutal domaine frontière d'un paraboloïde complet u€[0, ce] 

Application aux problèmes aux limites intérieur/extérieur de Dirichlet/Neumann inhomogène en v, homogène 
en u sans dépendance azimutal domaine frontière d'un paraboloïde tronqué u€ [u0, ce] 

Application aux problèmes aux limites intérieur de Dirichlet/Neumann, homogène en v, inhomogène en u sans 
dépendance azimutal domaine frontière d'un paraboloïde tronqué u£€[u0,c], condition inhomogène radiale, 
condition homogène angulaire 


Solution de l'équation de la chaleur plus généralement de l'équation de diffusion (régime transitoire) 


L'équation de la chaleur, séparation des variables d'espace et de la variable temporelle (à changer dans le texte) 
Configuration n°1 - Les conditions aux limites du domaine sont indépendantes du temps 
Cas N°1 : Les conditions aux limites sont homogènes 
Cas N°2 : Les conditions aux limites sont inhomogènes 
Configuration n°2 - Les conditions aux limites du domaine dépendent du temps et de l'espace : Une version du 
Théorème ou principe de Duhamel appliqué aux éguations de diffusion homogène (sans terme source) (voir page 
37 du Livre de H.S.Carslaw et J.C. Jaeger, Conduction of Heat in Solids) : 
Cas n°1 : température prescrite (condition de Dirichlet) à la surface dépendante de l'espace et du temps, ou flux 
prescrit (condition de Neumann) ou condition mixte de Robin 
Cas n°2 : (condition de Dirichlet) à la surface, ou flux prescrit (condition de Neumann) ou condition mixte de 
Robin, dépendante uniquement du temps, simplification du problème 


Configuration de domaine rectangulaire (cartésien) 
Exemple divers 


Configuration de domaine cylindrique : 
Corps en disque circulaire 
Exemple divers 
Corps axial symétrique (rz ) 
Exemple divers 


Configuration de domaine sphérique 
Exemple divers 
Cas de l'hypersphère : première partie, condition initiale et conditions aux limites de Dirichlet 
Cas de l'hypersphère : deuxième partie, condition initiale et conditions aux limites de Robin 


Configuration de domaine elliptique cylindrique 
Exemple divers 


Configuration de domaine parabolique cylindrique 
Exemple divers 


Solution du problème aux limites de l'équation de Poisson à deux et trois dimensions par la méthode de séparation 
des variables 


Méthode d'homogénéisation du problème par l'introduction d'une solution particulière 
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Les solutions élémentaires du Laplacien et les fonctions de Green d'une source ponctuelle placée dans un espace 
à N dimensions 


Bibliographie 
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Dans le vif du sujet : les équations de Laplace, Poisson et l'équation de la chaleur 


L'équation de Laplace en coordonnées cartésiennes s'écrit : 
O°T(x, y) F T(x, y) . 
ôx? ôy? 
L'équation de Poisson est proche de celle de Laplace et comporte un terme source : 
OT D), OT CRIER 
ox ôy 


0 








Ces deux équations désignent également les problèmes stationnaires de l'équation de la chaleur, 

qui dans sa forme temporelle peut se résumer à l'équation suivante comportant le terme de 

diffusion temporelle : 

T(x, y,t) y @T(x,y,t) 1 T(x, y,t) 
ôx? ôy? a ôt 





S(x, y) 


Dans notre cas nous commencerons la plupart du temps à considérer que les termes de sources et 
de diffusion sont des constantes. 


Conditions aux limites et conditions initiales 


Les problèmes auxquels nous recherchons des solutions sont dans la catégorie des problèmes aux 
limites d'équation aux dérivées partielles linéaires avec des conditions initiales dans un domaine 
fermé du plan, suivantes : 


T(x, y,0)= n(x, y) 


a T(x,7,9+ p Eht 


mm 


= f(x, y) 


os 





n normale à OS 


Lorsque l'équation représente un régime stationnaire seule les conditions aux limites définissent le 

problème 
T(x, y,t 

a T(x, y,t)+ß aatal 


mm 


= f(x, y) 


os 





n normale à OS 
a =1,f =0 > conditions de Dirichlet Ty D| = f(x,y) 
OT (x, 1) 


a =0, B =1 > conditions de Neumann 
on m 


= f(x,y) 


&,ß quelconque => conditions de Robin 
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Conditions aux limites, condition initiale homogènes, inhomogènes 


Une condition aux limites homogène est une conditions aux limites dont les valeurs fixées de la 
solution ou de sa dérivée première, ou d'une combinaison linéaire des deux est nulle à la frontière 
du domaine d'étude, et inversement inhomogène. 





Homogène 
soit T(x, Vha 
soit P Au ———-) =0 
on aS 
T 
soit æ T(x, y,t)+ B RCE ————) =0 
on n 





n normale à OS 


Les conditions aux limites sont inhomogènes inversement lorsque l'une des trois conditions est 
non nulle sur tout ou partie de la frontière du domaine. Lorsque d'un domaine est quadrangulaire 
alors l'homogénéité des conditions s'apprécie dans chaque direction du système de coordonnées. 


Une condition initiale est homogène lorsque la valeur initiale à t=0 de la fonction est nulle sur le 
domaine, et inversement inhomogène. La plupart du temps les conditions initiales sont par nature 
inhomogènes. 


Principe de superposition de l'équation de Laplace (et également régime stationnaire de 


l'équation de la chaleur) 


Lorsque la frontière du domaine S est la réunion de parties disjointes, sous-ensemble de S, alors la 
solution du problème aux limites peut s'écrire : 


a@s= |Jəs; 


85,n8$ =Ø 


Tœ y)= LT y) 





a; T(x, y)+ B; Len D fu 
T.(x, y) solution du problème aS; ij 
a, T, y)+ 2, Lad y) 0 
OS; ,i#j 





Soit donc le problème aux limites peut se simplifier en découpant la frontière en zone spécifique 
où une condition aux limites s'appliquent, à l'exception de toutes les autres zones où la condition 
aux limites homogène s'applique. 
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Principe de linéarité des solutions du problème aux limites stationnaires 


Du fait de la linéarité de l'équation aux dérivées partielles, la solution du problème aux limites 








T, (x, y) solution du problème æ T(x, y)+ B Len y) = f(x,y) 
n las 
T, (x, y) solution du problème à T(x, y) + B——— L 7) = g(x, y) 
os 
OT(x, y) 


>aT,(x,y)+bT, (x,y) solution du problème a T(x,y)+B— | =a f(x,y)+b g(x,y) 


On las 


Les principes de linéarité et de superposition permettent de simplifier la plupart des problèmes 
aux limites afin de la exprimer sous forme de combinaisons linéaires de solutions génériques. 





Séparation des variables 


La méthode de séparation des variables ne s'applique que quand la géométrie du domaine S 
présente une ou plusieurs symétries permettant d'envisager la recherche d'une solution simplifiée 
sous la forme de produit de fonctions de chacune des coordonnées du système. Pour être très 
concret, on propose de rechercher les solutions de l'équation de Laplace sur un domaine 


quadrangulaire cartésien sous la forme de 7 (7) =U(xW O). La séparabilité de l'équation est le 
fait d'aboutir à un système d'équations différentielles ordinaires couplées par des valeurs fixes. 


Il est possible de montrer que la séparabilité de l'équation de Laplace n'est possible que dans un 
nombre restreint de système de coordonnées curvilignes dit orthogonaux pour des domaines 
"quadrangulaires" de ces systèmes de coordonnées pourvu que ces même systèmes de 
coordonnées respectent une certaine symétrie des coordonnées induites par la forme de 
l'équation aux dérivées partielles. Toujours est-il que pour les systèmes de coordonnées, cartésien, 
polaire, cylindrique et une bonne dizaine d'autres la séparabilité est possible et permet d'envisager 
de trouver une solution aux problèmes posés. 


Volontairement nous employons u,v comme système de coordonnées pouvant désigner divers 
systèmes comme le cartésien x,y , le cylindrique r,ô, ou rz, ou z, ou le sphérique r, à, 

T(u,v,t) =U(u)(v)F (0) 

La méthode consiste alors à injecter cette solution dans l'équation aux dérivées partielles qui se 
révèlent être séparables c'est à dire réductible à un système de plusieurs équations différentielles 
ordinaires dépendant de paramètres particuliers, qui selon les types de conditions aux limites sont 
en nombre infini. Chaque équation différentielle est soumise aux conditions aux limites sur ses 
propres paramètres (u,v et t). Ce sont les conditions aux limites homogènes qui conduisent à une 
infinité de solutions qui se révèlent être les valeurs propres d'un certain opérateur différentiel (on 
se réfère en général à ce qu'on appelle l'étude spectrale de l'opérateur différentiel). 


La solution doit alors être recherchée comme une combinaison linéaire des solutions à valeur 
propre particulière. Et ce sont les conditions aux limites inhomogènes qui vont permettre de 
calculer les coefficients de la combinaison linéaire. 
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Problème de Sturm-Liouville régulier (et auto-adjoint) 


Parmi les équations différentielles que l'on rencontre le plus souvent, se trouve celle dîtes de 
Sturm-Liouville. 


Le Lever 0 6 Le TE ) fran) =0 


ZE [zz] p(z)>0 w(z)>0 








CL. a, Li +BU(z) =0 
ZO puo =0 
dx 





X2 

On va supposer que pour une telle équation différentielle les fonctions p(z) ont les propriétés 
suivantes sur l'intervalle [z:,z] : p(z) > O et w(z) > 0. Un tel système est communément appelé un 
système de Sturm-Liouville régulier pour des conditions aux limites mixtes homogènes. En voici les 
principaux résultats qui sont en général donnés sous forme de théorème : 

- on introduit un opérateur différentiel de Sturm-Liouville sous la forme : 


L(U(z))= l-le ROUE ot (pQU'(2))+sQU(2)) 


- l'opérateur différentiel de Sturm-Liouville est auto-adjoint sur l'espace de Hilbert réel avec la 


z 2 
f dz w(z) L(U(2))V (2) = f dz w(z)U(2) L(V (2)) 

mesure w(z) :4 z1 , ce qui a pour conséquence que les 

valeurs propres de l'opérateur différentiel de Sturm-Liouville sont toutes réelles, et encore bien 

d'autres propriétés (en nombre infini, oscillation des solutions, ...) 


- toutes les valeurs propres de l'opérateur différentiel de Sturm-Liouville L(U(z))=4U(2) sont 
réelles et non dégénérées (à une seule valeur propre correspond une seule fonction propre) et en 
nombre infini, rangé par ordre croissant avec pour limite l'infini : 

A <A, <... <A, <.. lim À, = +00 


n—> +0 n 


- Signe des valeurs propres : chaque valeur propre satisfait à l'égalité du quotient de Rayleigh : 


fear, oe oroe o) Eor or oF -+ febo oroe) 
À, =. z Z2 T Z2 5 
fa va, 2} fa wole, ©} 


Zi z 





PEDD, GP, E) -PED 0, G) ebole O ole, o) 
À, = Za + 
faz wo, DF 





Cette égalité permet d'assurer dans certains cas le signe des valeurs propres. 
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Une intégration par partie du numérateur de ce quotient donne les cas de valeurs positives des 
valeurs propres : 


hor 06, oF + febo, of roe, oF) 
4, = - 


n 





fa wo, DF 


jebol, oO oe oF) 
Cas 1 Dirichlet homogèneen z, etz, > D, (z)=®,(z,)=0= 1, =% >0 


fæ wolo, DF 


Zi 





ja (bolo, F + s(z)(© A (2}) 


Cas 2 Neumann homogèneen z, etz, = D, '(z)=®,'(z,)=0= 21, =< >0 


fa wolo, DF 





Cas 3 Mixte homogène en z, et z, Va,,@,,/,,B, 


CAN "()+89, (z)=0— D, "(Z)= Lo, (z,) aD '(z2)+ P-P, (z,)=0=> D, '(z,)= -Lo, (z,) 
pe) lo, GE) -rmo + fe (ae, +5, @)) 
l= 2 1 zi 


n 











fa wolo, 2) 
ou bien | 
COPA) EC UN CNE feat, a) +50, @}) 
je 2 — z 
faz walo, oF 
ou bien | 
pe) (6, GA} - pe) Žo, GŸ + Je. o toe, oF) 
2, — 2 a z 
fé we, OF 
ou bien | 
pe) loe- pe On EDP + feat, ooe, o }) 
des 2 1 z 





n 


fa walo, OF 


SE a aa aa ae a 
LE 1 1 2 
Cas 4 Neumann et / ou Dirichlet homogène en z, et /ouz, 


Sia,=0etB=0—2,>0 , Sia =0et B,=0— À, >0 
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- à chaque valeur propre est associée une fonction propre unique D,E) sur [z:,2] (à une 


D,(2) 


constante multiplicative près). Le système de fonctions propres est un système orthogonal 


sur l'intervalle [Z1,22] pour le produit scalaire de mesure w(z) soit 


dz w(z)®,(2)D,(2)=0 sin#m | | | 

ži . De plus il est complet, c'est à dire que toute fonction 
intégrable (avec le poids w(z)), continue et dérivable par morceau, peut se représenter par une 
combinaison linéaire infini des fonctions propres du problème de Sturm-Liouville. Cette dernière 
série fini converge sur tout l'intervalle : 


SOE CE, c= [7 a wL |? d wF 


-La norme des fonctions est calculée pour chaque problème de Sturm-Liouville sur le produit 
2 X2 2 
P, = |? d wP, (2) 
scalaire de mesure (2) : | Í. 
- le théorème d'Abel stipule que pour deux solutions (indépendantes) du problème de Sturm- 
Liouville, le Wronskien a la forme (ce qui garanti leur indépendance linéaire): 
U (z), V (2) 
(PU'E) s(z)+4w2) U) =0 (GA) 4-5) +22) (2) = 0 
PGYU'(2) (2)-U(2W'(2))= Cste W(U(2),V(2))=U'(2 (2)-U(2W' (2) 





Cste 
W(U(2),(2))= 
p(z) 
- Les fonctions propres ont un comportement dit "d'oscillation". Les théorèmes "d'oscillation" 


®, (2) 


stipulent que pour une fonction propre 
D, (2) 


, entre deux zéros successifs de la fonction propre 
: ; ; ` ; ; P, (zZ) avec A4 >A 
, il y a nécessairement au moins un zéro de la fonction propre 2 
; D R p, (z PRET ; | 
strictement entre les deux zéros successifs de a ( }. En général il y a plus de zéros de la fonction 


DANSE FANS AER que de la fonction Da GX soit à mesure de l'augmentation des valeurs 
propres. On dit que les fonctions propres "oscillent" de plus en plus vite à mesure de 
l'augmentation des valeurs propres. 

- on peut estimer la distance minimale entre deux zéros successifs dans le cas où p(z) ne s'annule 
pas dans l'intervalle [z:,z] : 


d rin (zéros successifs de D, (z))= A M, = malt s(z)+4,w(2))Sin?° (0) + 


À 


2 
Ce) ze [zz 19 € bar] 
p(z 


n 


- Formule asymptotique des valeurs propres : les valeurs propres sont espacés régulièrement 
asymptotiquement à l'infini se la loi : 


_2 
-rje fs PE ) lorsque n — œ 
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Problème de Sturm-Liouville régulier, auto-adijoint avec des valeurs propres tendant vers -°° 


Les raisonnements précédent restent inchangés dans le cas où les fonctions p(z) sur l'intervalle 
[zz] : p(z) < O et w(z) < 0. Le système de Sturm-Liouville reste régulier par évidence et en 
changeant tous les signes des expressions il vient que le système de valeurs propres est de signe 
inverse et tend donc vers --c. Attention il ne s'agit pas du cas où la fonction p(z) change de signe 
sur l'intervalle [z:,z,] pour lequel les conséquences de l'existence et de la répartition des valeurs 
propres sont légèrement différente. 


Exemple du cas des fonctions de Legendre dans le système sphéroïdal allongée. Le problème de 
Sturm-Liouville s'écrit alors : 


ARO e) + vs) =0 


zez,,z,| {z,,z,}e[l,oo| p(z)=1-2 s(z)=0 w(z)=1 


O O ERREUR) 
> l ) dz’ : dz ($+ vo fi ) dz | ($+ vo ni 


p&)<0 ze [zz] 








Notation A= {ee =v(v +1) = v= -2+ iT 


En inversant le signe de la fonction p(z), il vient donc que les valeurs propres tendent vers --c, dans 
ces conditions au vue de l'expression des valeurs propres, il vient alors des valeurs propres v 
intervenant dans les fonctions de Legendre sont imaginaires de la forme v=-1/2+it. Cela aboutit à 
l'utilisation des fonctions coniques de Mehler comme fonction propres (de première comme de 
seconde espèce) et de l'existence de valeurs propres de spectre discontinu lorsque l'intervalle du 
problème de Sturm-Liouville reste fini et que la fonction p(z) ne s'y annule pas. 


Problème de Sturm-Liouville singulier et auto-adjoint 


Les problèmes de Sturm-Liouville sont qualifiés de singulier dans les deux cas suivants : 
- la fonction p(z) s'annule sur l'une des bornes de l'intervalle /z:,z:] 
- l'intervalle [z:,z,] devient infini 


Dans le premier cas comme la fonction p(z) s'annule sur une des bornes de /z,z] et dans ce cas 
sur la borne en question, il n'est plus besoin de conditions aux limites fixées sur cette borne. On 
montre facilement que l'opérateur de Sturm-Liouville reste auto-adjoint. Dans le deuxième cas on 
suppose que l'expression Vw(z) f(z) tend vers 0 lorsque |z|->*°. Même si nous ne le montrons pas 
ici, les résultats sur l'existence, la réalité des valeurs propres et la complétude et l'orthogonalité des 
fonctions propres restent valides dans ces conditions. 
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a des fonctions de Legendre est typique de cas sur l'intervalle [-1,1] : 
dU 
Z» AE vo), awo -s@IU() = 0 
z e |-1,1] a =1-2? s(z)=0 w(z)=1 
2 
(re 2) Do), 0 pute 
=> dz dz 
p)>0 zel-1i] pŒD=p0)=0 
U(z)=P(z) v=n=0.12,.,n 
L'exemple des fonctions de Hermite sur l'intervalle [-s,1], solution du problème de Sturm- 


Liouville suivant : 
ZE |- 00 ,œ| p(z)= e7 s(z)=0 w(z)= e” A=2n 


AL DCE ve) +[1wG@)-sQU( = 0 








g? CUG) da zae” UO Fin UEO 
dz dz 
2 
se Z ar AD enO 
dz dz 
p(z)>0 ze |- %0,0] Lim p(z)=0 





Dans certains l'existence des valeurs propres est acquise, mais leur spectre devient continu (cas 
des problèmes aux limites sur des hyperboloïdes de révolution à une ou deux nappes. Dans ce cas 
l'emploi de transformation intégrale spécifique dont la forme est issue de la séparation des 
variables devient un outil précieux (voir transformation de Mehler-Fock). 
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Calcul des normes des fonctions propres du problème de Sturm-Liouville 


Lors de la démonstration de l'orthogonalité des fonctions propres à valeurs propres différentes, on 


parvient aux formules suivantes, si l'on prend deux fonctions propres : 





a o, (z 0® , (z f 
D, (2), (2) > (4 -4 )| dz wW2)®, C)9, (2) = Fe (22 Do D, (2)- 29, ©) 


Zi 
Passage à la limite et dépendance de la solution en À 


ôP, (2) 
=A +ô PD anD, GRR g 


R a(z 0®, (z k 
(a -4 )f dz w29, (2, (= pe (ee Do, (2- Le, @)] 


Z 
1 1 


3 oP, (z) 
= —ô} | dz w(2)®, ofo, (z)+ 64 mo) 


Z2 


op à (z) 


fo, (z)+64 ) 
A DG 0® , (z) Lo, GES 0®, 2) 
i ôz i oA 


Oz 





=| p(z) 


Zi 


8D, (z) _ 8P, (2) ôP, (2) 





<> a | w(z)®, (ZY = (61) Y [ae W(2)D, ()—"— Pae 2 z a refo, (z) re 


Zi Zi 





R 0®, (z) 0®, (z 8’, (z g 
ef snove ~ pef a aO te | +O(62) 


Zi 
1 


Il vient par passage à la limite 6À->0, la formule : 


wo, (z}? -| pe (2 = ôP (z) D, (z are) P, zao)| 





oÀ OO 


Zi Zi 


Exemple : cas des fonctions de Legendre : 


e problème de Sturm-Liouville s'écrit : 
dU E 
pO — A +A we) -sQIU() = 0 


PU(2) 


= +AU(z)=0 





p@)=1-2 s(2)=0 wz)=1= (1-2?) -r 


Notation À = p(p +1) = då = (2p + l)dp 


La valeur de la norme devient donc : 





® (2) 
(2p+D\ © Ôp z 


n020, 0> fau, EA a 


Zi 


ô® (z) D (z) zeo) 


Oz 


oÀ 


Ji 


1 
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Exemple : cas des fonctions de Bessel : 
Le problème de Sturm-Liouville s'écrit : 


a (z TE 2 }ave-stee =0 








ORO 2 2 OO O, 


Solution D'(2)+ J, WA z) v=4V2- (26 D!) 


La valeur de la norme devient donc (en simplifiant la notation) : 
A =v° >d =2vdv ©; (z) ®,(z) 


jo aee sate) 





Zi 


D'autre part les fonctions propres ont une propriété supplémentaire, à savoir 








sgela SO y à O eri 
V z 
v ôP, (z) = Z oP, (2) 
Ov Oz 
2 2 
Il vient 0P, (2) = á 3D, (z) => Q D, (z) = © Le + 1 0®, (z) 
v v & OvOz v © v © 


Porté dans l'expression de la norme cela devient : 


D, (= | dz w29, C} Jė (22 Zoe f GMO 2) 








Oz êz’ v Oz 


Zi 


_|_z ( (00) 7 = eo ee L 
T222) P A en || 


@’ D 0, p(z) 0®, (2) 
êz? Oz Ôz 


1 





1 





+(2w(2)-s(2)}D, (2) = 0 


Equation de Sturm — Liouville => p(z) 


0°®, (z) _ 1 (22 0®, (z) 
O7? o p(z) Ôz 


2| z (RO so 080.0 20 7 PO 4, nl 
læ, Lt s | DD +0, EP PE (a D |] 





+(4 w2) -s(2))®, Œ) 








Avec les fonctions de Bessel cela donne : 
2 


p(z)=z s(z)= _ w(z) = Zz 


2 


.e_1lf1(@æ@) vf | 
Eo | = EE = Jeo) í 5) 


Pour n=0 >| o| = i = 2) n eor) 
Z 


1 
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Exemple : cas des fonctions sinusoïdales : 
Le problème de Sturm-Liouville s'écrit : 


a (z TE D hane-stee =0 


d’U(2) 
dz? 











p(a)=1 s(2=0 w(z)=1= +AU(2)=-08 to +wU(z) = 0 


Cc 
Solution ®,(z)x Waz) v= >09, 90,2) 
in 


La valeur de la norme devient donc: 
A =v? = dÀ=2vdv 


fawo, > |1 ôP, (z) 0®, (2) MRE ,2 2,0) D, (z) 
p 2v oz Ov Ovoz J|, 
D'autre part les fonctions propres ont la même propriété que les fonctions de Bessel à savoir 
2 2 
y 0®, (z) 2e 0®, (z) = 0®, (z) iZ 0®, (z) á oP, (2) Z Ô P R 1 6®, (2) 
Ov z ôv v &z Ôv ©z v & v &z 
Porté dans l'expression de la norme cela devient : 


®, (2) > fd woa, e? IE (22 Z D,E f 0°? AOR 2) 




















Oz êz’ v Oz 


Zi 


-| (20) PP AO EUR AO) zno) 
2v Oz Oz 


CR) v@), p(z) 0®, (z) 
Oz? Oz Oz 


1 








+(2w(z)-s(2)}D, (2) = 0 


Equation de Sturm — Liouville => p(z) 


0°, (2) _ 1 (22 0®, (z) 
ð? p(z) Oz 


21 (0, Oo DO, OP) 2 PQ), A8, A j 
eol Lt T | Pre PRESS 2 0) FE ) 





+Q TOOLAG) 








Avec les fonctions sinusoïdales cela donne : 
p(z)=1 s(z)=0 w(z)=1 


= lo, @f =2 er. (2 ay) = voe 





Oz 
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Finitude des solutions 


Les solutions d'un problème aux limites ne sont pas qu'un simple problème formel mathématique 
déconnecté d'autres réalités, comme par exemple la modélisation d'un problème physique. Toute 
modélisation physique implique le plus souvent que les valeurs admises d'un problème purement 
mathématique soit finies. L'admissibilité de solutions divergentes en un quelconque point du 
domaine d'étude Q n'aurait que peu de sens. A contrario la finitude des solutions est une bonne 
hypothèse de base pour que toute solution mathématique puisse au moins représenter un 
phénomène physique. 


Singularité du système de coordonnées et régularité du gradient 


Cette notion est en quelque sorte une extension du principe de finitude des solutions de l'équation 
de Laplace. Il arrive que pour modéliser un problème physique , il soit plus commode d'introduire 
en sus des coordonnées cartésiennes d'autres systèmes de coordonnées, notamment ceux dits 
orthogonaux, afin d'appréhender des géométries particulières (cylindriques, sphériques, 
elliptiques, paraboliques, sphéroïdales, ...). Ces nouveaux systèmes sont décrits le plus souvent par 
des formules de correspondance avec le système cartésien. Dans bien des cas les formules de 
correspondance, et les grandeurs propres à décrire la nouvelle géométrie, introduisent des 
singularités qui ne sont pas présentes en géométrie cartésienne. Les singularités peuvent par 
exemple se traduire par l'annulation du Jacobien dans une certaine partie du domaine d'étude QO, 
des nouvelles coordonnées orthogonales. Dans le même temps cela s'accompagne d'une 
singularité apparente du gradient sur cette même partie. Les singularités du nouveau système de 
coordonnées peuvent également consister en une discontinuité dans leur valeur ou leur dérivée 
première lors d'un passage continu dans l'espace. 


Étant donné que cela n'est induit que par un effet de changement de variables, il n'y a aucune 
raison que la solution présente in fine la même singularité, une fois de retour dans le système 
cartésien. Ily a donc la nécessité que les solutions admissibles garantissent la régularité du 
gradient dans les sous-domaine où le Jacobien de la transformation du système de coordonnées 
s'annule. 


I n'est en général pas nécessaire de mentionner cette contrainte supplémentaire, car dans certains 
cas la condition de finitude suffit pour assurer la condition de régularité du gradient (exemple 
coordonnées cylindriques ou sphériques). C'est pourquoi on a du mal à comprendre au premier 
abord pourquoi la condition du gradient est demandée dans des systèmes comme le système 
parabolique cylindrique, ou le système sphéroïdale allongé ou aplati (prolate spheroïdal ou oblate 
spheroïdal en anglais). L'explication tient au fait que la condition de finitude garantit la régularité 
du gradient. Mais ce n'est pas toujours le cas, par exemple la régularité du gradient en système 
elliptique cylindrique est synonyme de continuité des solutions et de leur dérivée première au 
passage de la ligne dîtes inter-focale. 


Qu'à cela ne tienne, il est toujours bon de présenter les singularités du système de coordonnées et 
de s'assurer du respect de la condition de régularité du gradient. Nous tenterons à chaque fois 
d'évoquer le problème et de comprendre l'articulation avec le problème de la finitude de la 
solution. 


